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Abstract 

Based on Minkowski's work on critical lattices of 3-dinicnsional 
convex bodies we present an efficient algorithm for computing the 
density of a densest lattice packing of an arbitrary 3-polytope. As 
an application we calculate densest lattice packings of all regular and 
Archimedean polytopes. 



1 Introduction 

Throughout this paper, M*^ denotes the d-dimensional Euclidean space with 
origin 0, Euclidean norm || • ||, inner product (•, •) and unit sphere S'^~^. 
IC^ denotes the set of all convex bodies K d with nonempty interior 
int {K) and /Cg denotes the subset of /C^ consisting of all bodies which are 
centrally symmetric with respect to the origin. For a set M C M'^ we denote 
by vol(M) its volume with respect to its aflfine hull aff (M). Furthermore, 
conv(M), lin (M) denotes the convex hull, linear hull of M, respectively. 
The boundary of K S /C'' is denoted by bd {K). 

By a lattice KcW^ with basis S = {&^ . . . , b'^], where , . . . ,h'^ e 
are linearly independent, we understand the set 



d 



The determinant det A of A is the volume of the parallelepiped spanned by 
b^,... ,b'^, i.e., det A = | deti?|. A lattice A is called a packing lattice for 
K € /C^ if x + K and y-\-K do not overlap for y G A, x 7^ y. There always 
exists a packing lattice A*{K) with minimal determinant. Such a lattice is 
called a densest packing lattice and the quantity 

d*{K) =vol{K)/detA*{K) 

is called the lattice packing density of K or the density of a densest lattice 
packing of K. There is a large amount of literature on packings. Some 



*The second named author acknowledges the hospitality of the International Erwin 
Schrodinger Institute for Mathematical Physics in Vienna, where a main part of his con- 
tribution to this work has been completed. 



1 



books, which give a good description of the background of our work, are 
[|EGH89| , [lGL87|| , [|Rog64|] , iZon96|] as well as the Diploma the sis [lHau9g|] . 
For a more general survey on the theory of packings we refer to [FK9S] and 
the references within. 

The lattice packing problem for a general body in IC^ is very hard. In 
fact, for d > A the only exact results are on space fillers (cf. [tVIcMSq] , 
[|Ven54| ) for which 5*{K) = 1 and on the unit bah B'^ where 5*{B'^) is 
known for d < 8 (see e.g. ]CS98| ], |Zon99[ ) . In contrast for a fixed body 
K ^ K? there are several techniques to solve the problem (see e.g. |GL87], 
pp. 241) and there exists also an algorithm, due to Mount and Silverman 
[ MS90| , that determines 6*{P) for a centrally symmetric n-gon in time 0{n). 

However, already in 3-space the situation is rather more complicated. 
Apart from the 3-dimensional space-fillers (for a classification see [GL87], pp. 
164) and from cylinders based on a convex disk, in which case it can be shown 
that the problem is equivalent to the determination of the lattice packing 
density of the convex disk (cf. [Rog64], p. 13), densest lattice packings are 
only known for the bodies in M^, which are listed in the following table: 



Body 

regular octahedron 

truncated cubes; for < A < 3 : 

{x G : \xi\ < 1, |xi -I- 2:2 + 2;3| < A} 

|x e : (xf + xlY'"^ + Ixsl < l| 

frustrum of a sphere; for < A < 1: 
[x : xl + xl + xl < 1, IX3I < A} 

tetrahedron and cubeoctahedron 



Author 

Minkowski [|Min04|], 1904 



Whitworth iWhi4g|] , 1948 



Whitworth jWhiSlI] , 1951 
Chalk [phaSqi, 1950 



Hoylman [|Hoy70|| , 1970 



It is worth to mention that the family of frustrums of a sphere includes 
the 3-ball as a limiting case, for which the packing density was determined 
already by Gauss []Gau31 |. 

All the computations of 5*{K) for the bodies in the table above may 
be regarded as an application of a general method developed by Minkowski 



[Min04| which characterizes densest packing lattices of a 3-dimensional con- 
vex body by certain properties (see also [|GL87 |, pp. 340). However this 
method was considered as rather impractical and in 1964 Rogers [Rog64] 



wrote ^''Despite considerable theoretical advances in the Geometry of Num- 
bers since Minkowski's time, the problem of determining the value of 5*{K) 
for a given convex 3-dimensional body K remains a formidable tasJi\ In- 
deed the only change in the list of known densest lattice packings since the 
publication of Rogers' book is the addition of the tetrahedron and cubeoc- 
tahedron. And even in 1990 Gruber mentioned the determination of 5*{K) 
of a 3-dimensional convex body as one of the important open problems in 



Geometry of Numbers | Gru9C ] . 



2 



Here we use Minkowski's work as a starting point for the construction 
of a practicable algorithm to compute the packing density of an arbitrary 
polytope in M^. We proceed as follows. In section 2 we adapt Minkowski's 
work to our purposes. Moreover we discuss in this section the principal 
applicability of his work in higher dimensions. 

Having adapted Minkowski's method we still have to solve two prob- 
lems to obtain a practicable algorithm. First we have to determine the 
local minima of a polynomial of degree three in three variables. Somewhat 
surprisingly there appears to be no general purpose algorithm of numerical 
analysis which suits our needs. Thus we develop an ad hoc method for our 
problem. This will be done in section 4. Before doing this we deal in sec- 
tion 3 with a more geometric problem. It turns out that by the methods 
of section 2 we have essentially to look at every choice of 7 facets of the 
polytope. This leaves us with more or less (") cases for a polytope with n 
facets and in every step we have to solve a nonlinear minimization problem. 
This limits the applicability of the algorithm to polytopes with few facets. 
Thus we develop in section 3 some methods to reduce the number of cases. 
While we do not give an exact worst case analysis it should be possible to 
reduce the number of cases to possibly as few as 0(n'^/^) in typical cases 
by showing that optimal packings must lead to feasible points of some re- 
lated simple linear optimization problems. This will be done in section 3. 
It should however be mentioned that in our program we do not fully exploit 
this reduction as an implementation would become rather complex. As an 
application of our work we present in section 5 optimal lattice packings for 
all regular and Archimedean polytopes. 



2 Necessary conditions for optimal lattices 

In this section we state without proof Minkowski's results. In fact, we have 



summarized the relevant facts in Theorem 2.1. We have included some more 
results, as we discuss at the end of the section, whether the algorithm could 
be extended to higher dimensions. Further we have changed and modernized 
the notation. We remark that a good part of Minkowskis theory is exposed 
in [ |GL87 | and |Zon9(;], though in both books Theorem is stated in a 



slightly weaker form. 

For Ki e /C^, Xi £ R, i = 1,2, we denote by XiKi + X2K2 the set 
{Xix^ + A2X^ : £ Ki, G K2}. Minkowski observed, that 

A is a packing lattice of K ^ h. is a packing lattice of — K). (2.1) 

Since the difference body ^{K — K) belongs to the class /Cq, in the following 
we assume that all bodies are centrally symmetric. 

A lattice A is called admissible for K E /Cq, if int X n A = {0}. It is well 
known and easy to see that A is admissible if and only if 2A is a packing 
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lattice. The value 

A{K) = minjdet A : A admissible for K} 

is called the critical determinant of K and an admissible lattice A satisfying 
A{K) = detA is called a critical lattice. Thus 6*{K) = vol{K)/{2'^A{K)) 
and the problems of constructing a densest packing lattice and a critical 
lattice arc equivalent. 

While we naturally do not know a basis of a critical lattice beforehand, 
we shall see that we have a great amount of information on the behaviour 
of certain lattice points with fixed coordinates with respect to such a basis. 
For a basis B = {b^,. . . , b'^} of and a point x e we denote by xb = 
{xi,... , Xd)^ the vector given by xb = Y.i=i ^i^^- 

The construction of critical lattices is based on the connection of lat- 
tice crosspolytopes and primitive vectors. As usual, a set of lattice vectors 
b^, . . . ,b'^ of a lattice A C M*^ is called primitive iff this set can be extended 
to a basis of the lattice A. Let A C M'^ be a lattice and let 6^, . . . , ft'^' G A be 
linearly independent. The crosspolytope C = conv {±b\...,±b''} is called 
a lattice crosspolytope iff A PI int (C) = {0}. If we even have A n bd (C) = 
{±6^, . . . , ±6^^} then it is called a free lattice crosspolytope. Clearly, the con- 
vex hull of every primitive set b^, . . . forms a free lattice crosspolytope 
and every free lattice crosspolytope is a lattice crosspolytope, while the con- 
verse is not true. A complete characterization of (free) lattice crosspolytopes 
of dimension up to three is given in 

Lemma 2.1 (Minkowski, 1904). Let A be a lattice in M^. 

• For k = 1 the vertices of any lattice crosspolytope are primitive. 

• For k = 2 the vertices of any free lattice crosspolytope are primitive, 
while for the non-free lattice crosspolytopes C there exists a basis B of 
the lattice such that C = conv {(1, 0, 0)^, (1, 2, 0)b}- 

• For k = 3 there are two types of free lattice crosspolytopes. One has 
primitive vertices (crosspolytope of the first type), the other has the 
vertices (1,0, 0)^, (0, l,0)s, (1, 1,2)b for a basis B (crosspolytope of 
the second type). For every non-free lattice crosspolytope C there 
exists a basis B such that C = conv {(1, 0, 0)b, (0, 1, 0)b,p}, where p is 
element of the set {(0, 1, 2)^, (1, 1, 2)^, (1, 2, 2)s, (1, 1, 3)s, (1, 2, 3)s, 
(2,2,3)b, (1,2,4)b, (2,3,4)b}. 

Using free lattice crosspolytopes we can identify critical lattices for every 
K e )Cq. To this end we use the following abbreviation: For a basis B = 
{b\b'^,b^} of M3 let 
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Uh = {{1,0,0)3, (0,1, 0)b, (0,0,1)b, (0,1, -1)b,(-1,0,1)b, (1,-1, 0)b}, 
Z^| = {(1,0,0)b, (0,1, 0)b, (0,0, 1)b, (0,1, 1)b, (1,0, 1)b, (1,1,0)b}, 
Z^| = {(1,0,0)b, (0,1, 0)b, (0,0,1)b, (0,1, 1)b, (1,0,1)b, (1,1,0)b, 
(1,1,1)b}. 

(2.2) 



Lemma 2.2 (Minkowski, 1904). Let K £ IC^. Then th ere exists a criti- 
cal lattice A with basis B such that one of the following cases holds: 

1. U]^ C hA{K) and K contains no lattice crosspolytope of the second 
type. 

2. C bd (K) and (1, 1, l)ij ^ K. 

3. Ufj C hd{K). 

While the distinction between cases 2. and 3. may look artificial, we shall 
see that it leads to rather different cases in the actual computation of a 
critical lattice. For any lattice with basis B which satisfies a condition of 
the previous lemma we can easily check its admissibility: 

Lemma 2.3 (Minkowski, 1904). Let G /Cq and A be a lattice with 
basis B. Then A is admissible, if one of the following conditions is satisfied 

1. Ul C bd(K) and (-1,1,1)b,(1,-1,1)b,(1,1,-1)b 0intK. 

2. Ul C hA{K) and (1,1,1)b 0intK. 

Proof. The second statement follows immediately from Minkowski's charac- 
terizations of lattice crosspolytopes (cf. [ ^on96| , Lemma 4.9, 2*), whereas 
item 1. is contained in a more implicitly way in his paper. However, from 
the work of Minkowski we know (cf. ||Zon96 |, Lemma 4.9, 1*) that in the 
first case, i.e., U]^ C bd(i^), the lattice A is admissible if 

al) (0, ±1, ±1)b i K, a2) (1, ±1, ±1)^ ^ K, a3) (2, ±1, ±1)b i K, 

(2.3) 

as well as all points arising from permutations of the coordinates of the above 
points are not contained in K. Therefore, in order to prove 1. we have to 
show that the points (0, 1, 1)b, (1, 1, 1)_b, (2, ±1, ±1)^ and the corresponding 
permutations are not contained in K. To this end let / : M'^ — > M>o be the 
distance function of K, i.e, f{x) = min{A : A G M>o and x £ \K}. Then 
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K = {x£R^ : fix) < 1} and since K £ ICq the function / describes a norm 
on M^. Since U}^ C bd (K) we find 

/((0, 1, 1)b) > /((0, 0, 2)b) - /((O, -1, 1)b) = 1, 
/((1,1,1)b) > /((3,0,0)b) - /((1,-1,0)b) -/((1,0,-1)b) = 1, 
/((2,1,-1)b) >/((2,0,0)b)-/((0,-1,1)b) = 1, 
/((2,-1,1)b) >/((2,0,0)b)-/((0,1,-1)b) = 1, 
/((2, -1, > /((2, -2, 0)b) - /((O, -1, 1)b)) = 1. 

Obviously, the same inequahties hold for the points given by all permutations 
of the coordinates of the points of the left hand side and thus 

/((2,1,1)b) > /((2,2,0)b) -/((0,1,-1)b) > 1. 

□ 



(2.5) 



We call the sets Z^^, j = 1,2,3, (cf. ( |2.2|) ) test sets of the first, second or 
third kind, respectively. 

Now suppose that B = {b^,b^,b'^} is a basis of a critical lattice A of -ftT 
and let Ub = {u^b^ ■ ■ ■ j'^b}^ A; = 6 or k = 7, be one the three test sets such 
that Ub C hd{K) (cf. Lemma 2^). Let Hi be any supporting hyperplane 
of K containing u^, i = 1, . . . ,k, and let 

Sh,,...,h, = {W£ m3x3 -.uIj^t e Hi, l<i<k} . (2.4) 

On this space we consider the function 

fHi,...,Hk ■ 'SHi,...,Hk K given by 

fH„...,H,{W) = \det{W)\. 

If we assume for a moment that Ub = A n bd (K) then 

S is a local minimum of /hj^...^j7j.(W), W G SHi,...,Hk- (2-6) 

Otherwise there exists a W G SHi,...,Hk ^ sufficiently small neighbour- 
hood of B such that | det(W^)| < det A and K n {WZ'^)\{0} C {u]y, . . . , 
u^f}. However, this set of vectors is contained in the supporting hyper- 
planes Hi, . . . , Hi; for any element W £ S{Hi, . . . , Hf^.). Thus the lattice 
WIj'^ is admissible and hence A can not be critical. 

In general the situation is much more complicated as Ub hd{K) may 
just be a proper subset of A n bd(i^). But by a close examination of all 
possible cases Minkowski found 

Theorem 2.1 (Minkowski, 1904). Let K £ }Cq. Then there exists a crit- 
ical lattice A of K with basis B such that bd (K) D A contains a test set 
^B ~ {^B' ■ ■ ■ ' '^b} f'^^ 0, 3 ^ {Ij 2, 3}, such that for any choice of support- 
ing hyperplanes Hi of K containing n^, 1 < i < k, one of the following 4 
cases holds: 
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I. j = 1 and ( p.6| ) holds. 
II. j = 2 and (U) holds. 
III. j = 3 and {^) holds. 



IV. J = 3 and i/iere are scalars Ai,A2, Ag > 0, A3 G M, swc/i i/iai /or the 
outer normal vectors v^,v'^,v^,v^ of the hyperplanes Hi, H2, H^, Hq, 
respectively, holds: a) Xiv^ + X2v'^ + X^v^ = Xqv^ andh) the hyperplane 
Hq with outer normal vector Xiv^ + X2v'^ containing is a supporting 
hyperplane of K and ( |2.6D holds with Hq replaced by H^ . 



At this point we should remark that in case I. Minkowski gave a some- 
what stronger condition, but our form appears to be more suitable for au- 
tomatic computation. 

Of course, given an arbitrary convex body K we do not know how to ex- 



ploit Theorem 2.1, but if we consider only polytopes then for the supporting 



hyperplanes Hi in Theorem 2.1 we may always choose the supporting hyper- 
planes of the facets of the polytope. As a polytope has only finitely many 
facets we obtain the following frame of an algorithm for the computation of 
a critical lattice of a polytope: 

Algorithm 2.1. Let P G /Cj]. 



For each of the cases I., II., III., IV. of Theorem 2.1 do 

— For every choice of k facets with supporting planes Hi , . . . , Hk of 
P (k = 6 in the first two cases and k = 7 in the latter ones) do 

51. Determine Shi,..., Hi:- 

52. Find the local minima M-Hi,...,Hi: 0/ /hi,...,//^ (c/. (|2.6|)). 

53. For each M G M.Hi,...,Hf, check whether M -lA is an admis- 



sible lattice, i.e., M satisfies the criterion 1. of Lemma "LZ 
in the first case and criterion 2. in the remaining cases. 

• Among all calculated admissible lattices find one with minimal deter- 
minant. The corresponding lattice is a critical lattice of P. 

It turns out that at this point we are left with two problems. First there 
appears to be no general purpose algorithm to find all local minima of a 
function like fHi,...,Hk ^ have to do in Step S2. Moreover, a priori we 
can not assume that the local minima of this function are isolated points. In 
general, they may form a manifold and we have the problem to parameterize 
such a manifold in order to carry out step S3. In section |^ we shall show 
how one can overcome these problems. Another problem is just the number 
of steps of the algorithm. In a straightforward implementation, we have to 
consider every choice of 6 and 7 facets of the polytope and hence we have 
about (") steps. Of course, this limits the algorithm to polytopes having 



7 



only few facets. Hence in order to get an efficient algorithm we have to 
reduce the number of steps. This can be done very effectively by some 
considerations given in section ^. 

We close the section with some remarks on the extension of the algorithm 
to higher dimensions. While Minkowski settles his work in 3-space, the 
ideas principally work in higher dimensions as well. In fact there has been 
an enumeration of lattice crosspolytopes in dimension 4 ( [|Bru44| , |Lau4f:], 
[[Wol54| ], ||Mor60|| , ||Ban65| ], ||Ban671] ) and dimension 5 ( [[Wes89| ]). Beside this 
classification we have to determine the number of different test sets. For the 
cardinality of a test set we have the general natural lower bound of d{d+l)/2 
given in [^wi53| ] . There is no obvious upper bound for their cardinality, but 
the results in dimension two and three suggest the upper bound of one half of 
the maximal number of lattice points contained in the boundary of a lattice 
point free strictly convex set. Due to a result of Minkowski this number 
is bounded by 2*^ — 1. Finally, we have to take into account the additional 
lattice points (not contained in the test sets) lying in the boundary, which 
are responsible for the split of test sets of the third kind into two separate 
cases in dimension 3. Thus even in dimension 4 there should be a large 
number of cases which have to be considered separately and there has been 
no attempt to give an enumeration of these cases. 

Further in each case we have to consider all possible choices of facets 
and even though we could apply the results of the next section it appears to 
be computationally difficult to choose at least 10 out of n facets. Moreover, 
even in dimension 4 we have to solve a minimization problem for polynomials 
of degree 4 with up to six variables. Again it appears to be not an easy task 
to find all local minima reliably and quickly. Thus without introduction of 
new ideas we have the impression that the algorithm is practically restricted 
to 3-space. 



3 Necessary conditions for test sets 



In the following let P € /Cq be a centrally symmetric polytope with n facets 
Fi and let Hi = aff (Fi), 1 < i < n. By we denote the halfspace bounded 
by Hi containing P and let H^^ = M.^\H^ U Hi. We always assume that we 
have a lattice description of the polytope, i.e., we know the face lattice of P 
specified by its Hasse diagram and the vertices and edges of P (cf. |Sei97| ]). 
In particular, for each facet Fj we have a list M{Fi) of its edge-neighbours, 
i.e., 

Af{Fi) = {Fj : dim(Fi n Fj) = 1}. 
We remark that such a lattice description can be computed from the descrip- 



tion P = D^^^H^ in time O(nlogn) |Cha93]. Regarding the combinatorics 



of polytopes we refer to the books | McSh71 ] and | Zie95|] 
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As pointed out in the last section one crucial point of Algorithm 2.1 
is the number of choices of facets (or hyperplanes Hi) which have to be 
considered for each case of the algorithm. In this section we show how one 
can reduce this number. However, since it turns out that the most time 



consuming step of Algorithm 2.1 is step S2., we are also looking for ways 
to reduce the number of executions of step S2. as well. An exact worst 
case analysis of the complexity of the resulting algorithm appears to be 
not completely straightforward. But we show that for some rather natural 
classes of polytopes with n facets we can eliminate "most" possible choices 
of facets in time O(n^) and we have to carry out only 0(n^/^) times the 
steps SI., S2. and S3. 

Of course, using the central symmetry of the polytope and the arbi- 
trariness of the order of the basis of a lattice we can reduce the number 
of possible choices of hyperplanes Hi in the first two cases to j2 (g) and to 
2 • choices in the remaining cases. This does not really help and so 

one could try to make further use of the symmetries of a given polytope, as 
Minkowski did in his study of the octahedron, where he managed to reduce 
the number of cases to 1 (!). Thus he did not need to carry out one step S2. 
However, for polytopes with little symmetry this would not be of great help 
and therefore we use a different approach. 

With respect to step S3, of algorithm |2.1] we are only interested in a 



selection of hyperplanes Hi, 1 < i < k, (say) such that 

Shu...,h, n{W e M^^^ : Uw Chd (P)} / 0, (3.1) 

where Uw is a test set of the first, second or third kind corresponding to 
the case we are studying. If Uw = {uw, ■ ■ ■ ,u^^} then ( p.lj ) just says that 
the point u\y should not only lie in the hyperplane Hi, but in the facet Fi. 
Hence (|3.1|) can be reformulated as the condition that 



(3.2) 



for all Fij G M{F'), 1 < i < A;} / 0. 



To check whether Shi,...,Hk is empty is just an instance of a feasibility prob- 
lem of Linear Programming and may be easily solved by any LP-solver (cf 
e.g. pch8^ ]). 

Since the vectors in a test set are linearly dependent the set Shi,...,Hk '^iil 
be empty for "most choices" of hyperplanes Hi. For instance: Let a = —1 
if we are dealing with case 1 and let o" = 1 otherwise. Then the vectors of a 
test set satisfy the relations 

If we have fixed a facet Fi with hyperplane //i = {x G : [a^yx = bi}, say, 
for the vector lijy, then by the first relation of (|3.3|) we get a{a^yu^^ < 0. 
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Otherwise the sum of these two vectors would be separated by Hi from the 
polytope, but the sum u^^ has to he in a facet. This trivial observation 
already reduces the possible choices for H2 to almost n/2 and obviously we 
can apply the same argumentation to the other hyperplanes. However, as 
we shall see at the end of the section, a detailed analysis of the linear depen- 
dencies will give a much better reduction for certain classes of polytopes. 

Of course, the determination of the emptiness of ( p.2| ) would reduce the 
number of executions of step S2, but we still have to consider all possible 
choices. Hence we have the problem to find a "fast" way to exclude "almost" 
all choices of facets (hyperplanes Hi) with Shi,...,Hk = ^- To this end we 
make use of (^) . Let 

g = {{Fi, Fj,Fi,) : {Fi + aFj) n ^ 0} . (3.4) 

Obviously, 

Sh,^ / =^ (Fi, , Fi, ,Fi,), {Fi, , Fi^ ,Fi,), {Fi, ,Fi^,Fi,)eG, 

and to test whether a tuple [Fi, Fj, Fk) belongs to G is just a feasibility 
problem of Linear Programming, namely: 

{F^,Fj,Fk) G g ^{iw\w^) G M3><2 : e H„ G H+ for Fi, G AA(F,), 

w'^ G Hj, G for F,, G AA(F,), 

+ aw"^ G Ffc, + aw^ G for F^, G AA(Ffc)} / 0. 

(3.5) 

The construction of the set Q by enumeration clearly involves the con- 
sideration of O(n^) possibilities. While in the numerical examples given 
in the last section the most time consuming part of the algorithm was the 
solution of (^.2|), the series of examples at the end of the section indicates 
that for large polytopes the determination of Q could be the hardest part. 
Therefore we show now how the geometry of the polytope can be used to 
do this effectively. To this end let 

Q[Fi) = {F, : [Fi + ctF,) n bd (P) / 0} , \<i<n. 

With this notation we have 

Lemma 3.1. Let Fi be a facet of P and let F\^ = \Jp.^g(^p.^ Fj . The set 
F^ is edge-connected, i.e., any two points x,y £ F-^ can be connected by a 
continuous path contained in F]^ without crossing a vertex of P. 

Proof. For a fixed point v G Fj let l{v) = bd (P) n o-(bd(P) - v) = 
bd(P)n(bd(P)-cJu). Then we have T(i;) = {y e Pr\{P - av) : y + ^iv, ^jl e 
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M, is a supporting line of PCi (P — av)} and thus I{v) is the shadow bound- 
ary of Pn (P — av) in direction v. Hence two points oiI{v) can be connected 
by a continuous path. Now let 

I(v)^ = U{Fj:F,nX{v)^9)}Fj. 

Then, by construction, the set I{v)'^ is edge-connected and so it is the union 

Once we have found one element of the set Q{Fi), the lemma says that we 
can determine the other elements by recursively checking the neighbours of 
the elements which were already found and since 

(Fi, Fj,Fk) £g^Fj£ g{F,), (3.6) 

we may use Lemma |3.1| in order to construct the set Q. In order to present 
an algorithm for computing Q as well as the sets Q{Fi) we need one more 
notation: For /,n € M^, Z < u, let B{l,u) = {x € : / < x < u} be the 
box parallel to the coordinate axes with lower vertex I and upper vertex u 
and for a facet Fi let B{P,u^) be the minimal box containing Fi, which can 
be computed from the coordinates of the vertices of Fi in time 0{n). A 
necessary condition for [Fi + crFj) n -F^ 7^ is given by 

{B{l\u')+aB{P ,u^)) n B{l'',u'') = 

B{1' + P,u' + u^)nB{l^,u^) a 
B{1' -u^,u' -P)nB{l^,u^) a 

which can easily be checked since 

B{1, u) n B{1, ii) / / < tt and [ < u. (3.8) 

Therefore, if we want to find for a given pair Fi,Fj all facets -F^ with 
{Fi, Fj, Fk) G G we just have to consider the facets corresponding to boxes 
B{l^,u^) satisfying ( [3.71 ). For polytopes with many facets "most" facets 
will be "far away" from Fi + Fj and thus ( tj.7| ) won't be fulfilled for "most" 
facets. 

The next lemma says how we can find one facet lying in a set Q{Fi) with 
the help of these boxes . 

Lemma 3.2. Let v{P) = max{#M{Fi) : I < i < n}, let f5{P) be the 
maximal number of boxes B{l^,u^) intersecting a box of the form B{1^ + 
u'^) + aB{P + u^), 1 < i,j < n, and letr]{P) be the maximal number of facets 
of P having a nonempty intersection with a fixed hyperplane containing the 
origin. Then a facet Fj lying in a set G[Fi) can be found in time 0{n + 
rj{P)log\n)+rj{P)(3{P)u{P)). 



= 1, (3.7) 

= -1, 
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Proof. Let v £ Fi and let H be any hyperplane containing v and the origin. 
Then there exists a facet Fj G Q{Fi) having a nonempty intersection with 
HnP. Using the edges of P we can easily determine all facets {Fj-^ , • • • , Fji}, 
ji ^ '>l{P)^ with this property. Since the polytope has 0{n) edges this can 
be carried out in time 0{n). For each j'j the boxes B{l^,u^), k < P{P), 
intersecting [B{1'', u*) + aB{P% u-'*)) can be determined in time 0(log^(n) + 
f3{P)) by well know methods from computational geometry about range 
searching (cf. [ [Aga97] ]). For each possible choice Fi,Fj.,Fk we use (3.5) 
to verify whether (Fi, Fj-, Fk) G G and thus Fj- S GiFi). Now each (3.5) 
is a feasibility problem of Linear Programming with 0(i^(P)) constraints 
in dimension 6 and this can be solved with 0(i^(P)) arithmetic operations 



(cf. gch8|], pp. 199). □ 



Using the last lemma we have the following algorithm for computing the 
set G and the sets G{Fi), 1 < i < n. 

Algorithm 3.1 (Computing G{Fi) and G)- 

Input: A polytope P G /Cg given by the supporting hyperplanes Hi, 1 < i < n, 
of its facets Fi, as well as a lattice description of P. 
Output: G and G{Fi), 1 <i < n. 

• Let [/ = and determine the boxes B{l^,u'^), 1 < i < n. 

• For each facet Fi do 

Gl. Find a facet Fj G G{Fi). Let Ui = % andN = {Fj}. 

G2. While N do 

G3. For Fj^ £ N determine all facets F^^ , . . . , such that 
{B{l\u^) + aB{V\u^^)) r\ B{l^i ,u^i) /0 (cf.^^^), (P)). 

G4. Use ( p.5D to determine all {Fi, Fj-^, Fk^) £ G, ki £ {ki, . . . , kp} , 
and add them to the set U. 

G5. Let Ui = UiU {FjJ and N = {N U M{Fj^))\Ui, if there exists a 
{Fi,Fj^,Fkt) e G, otherwise let N = N\Fj^. 

Lemma 3.3. Let j{P) = max{#^(Fj) : I < i < n}, and let v{P), r]{P), 
P{P) as in Lemma |3.2| . Then Algorithm \3.^ determines G and G{Fi), 1 < 
i<n,in time 0{n'^ + nlog^{n){rj{P) + -f{P)u{P)p{P)) + nP{P)u{P){r]{P) + 



Proof. By Lemma 3.1 the algorithm finds all {Fi, Fj, F^j G G and at the 
end we have U = G. Furthermore, at the end of each loop G2. the set Ui 
coincides with G{Fi). 

To find a first facet Fj G G{Fi) we need by Lemma p.2| at most 0{n + 
r]{P)log'^{n) + r]{P)l3{P)v{P)) operations. For the estimation of the steps 



G3. and G4. we can proceed as in the proof of Lemma 3.2. The boxes in step 



12 



G3. can be found in time 0(log^(n) + and in step G4. we have to solve 

at most (3{P) feasibility problems (cf. (|3.5| )) with at most iy{P) constraints in 
dimension 6 which can be done with 0(i/(P)) arithmetic operations. Finally, 
we observe that for each facet Fi the loop G2. is executed at most 7(P)z^(P) 
times. □ 

Remark. 

i) The bound on the running time given in the last lemma is useless for a 
worst-case analysis, because there exist polytopes such that each of the 
numbers r/(P), /3(P), 7(P), z^(P) is of order 0{n). In this case Lemma 
|3.3| would give a 0(n^)-algorithm and of course, one can determine 
the sets Q,Q{Fi) by a trivial 0(n^)-algorithm. Nevertheless we shall 
see in Theorem that this lemma gives a 0(v?) bound for a rather 
natural class of polytopes. 

ii) As we have Fj £ G{Fi) <^ Fi £ Q{Fj) we may use the facets belonging 
to a set Q{Fi) as starting facets in step Gl. It is not hard to see that 
we can determine all starting facets except the first one in this way. 
However, for simplification we do not exploit this fact. 



Altogether the previous observations lead to the following refinement of 



Algorithm 2.1 



Algorithm 3.2. 

Input: A polytope P £ ICq given by the supporting hyperplanes Hi, 1 < i < n, 
of its facets Fi, as well as a lattice description of P. 
Output: A densest packing lattice. 

• For each of the cases I., II., III., IV. of Theorem do 

II. Compute the sets G{Fi) and Q with Algorithm |3.1| . 

— For three facets Fi^,Fi^,Fi^ satisfying Fi^ G o,f^d Fi^ G 

g{Fi,)ng{Fi,) do 

— For every choice of facets Fi., 4 < i < k, with {Fi-^, Fi^, Fi^), 
{Fi„Fi,^Fi,),{Fi„Fi,,Fi,J do 

SO. IfSH,^,...,H,^^9id. O) do 

51. Determine Shi_^,...,Hi^- 

52. Find the local minima Mhi^,...,Hi^ of fHi^,...,Hi^- 

53. For each M A4h, h, check whether M ■ is an ad- 
missible lattice, i.e., M satisfies the criterion 1. of Lemma 



2.3 in the first case and criterion 2. in the remaining cases. 



Among all calculated admissible lattices find one with minimal deter- 
minant. The corresponding lattice is a critical lattice of P. 
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It seems to be a non trivial problem to give a "non trivial" worst case 
analysis of the algorithm for an arbitrary polytope. In the following we want 
to demonstrate, by a rather natural series of polytopes, the improvement of 
Algroithm (^]^) compared with the brute force method (see Algorithm ( |2.lD ) 
which involves the examination of Q{'n7) steps SO., SI., S2., and S3. 

For a facet F of a polytope P we denote by R{F) its circumradius and 
by r{F) its inradius with respect to its affine hull. For a real c > we say 
that the facets of P are of c-uniform shape, if 

min{r(F) : F facet of P} > c • max {i?(F) : F facet of P} . 

Theorem 3.1. Let {Pm}m€N be a series of polytopes Pm S ICq such that 
i) all Pm have facets of c-uniform shape for some fixed c and ii) {P™} 
converges to the 3- dimensional unit hall with respect to the Hausdorff metric. 
Let fm be the number of facets of the polytope Pm- Then the number of all 



choices of facets which will be examined by Algorithm (3^) in the steps 
SO., SI., S2. and S3, is of size 0{fm'^) and these possible choices of facets 
can be determined in time 0{fm) 



Proof. We project hd{Pm) by the radial projection p onto the unit sphere 
5^ and we carry out all calculations on S*^. It can easily be checked that our 
asymptotical estimates remain correct. The facets of Pm will be denoted by 
F/") 1 ^ i ^ fm, and in the following we shall denote by Cj certain positive 
constants. First we observe that the uniformity of our sequence implies that 
for the spherical diameter d{p{F^)) and the spherical area A(p(Fl")) of the 
facets holds 

ci/v^ < dipiFn) < C2/^/^^, cs/fm < A{p{Fn) < c^/fm- (3.9) 

Next we note that for a point x £ S'^ we have 

{y(.S^:y + ^^S^} = {yGS'': {x, y) = -1/2}. (3.10) 

Now let F^, F™ be two facets of Pm such that {F^ + cjF™) n bd {Pm) / 
and let x e F™. We find by (|]9|), (|I^) that 

F™ C {y G bd {Pm) ■■ -a/2 - cg/ < {p{y),pix)) < -a/2 + c^/^}. 

(3.11) 

As the spherical area of the set on the right is bounded by CQ/\/y^ we 
find by (^) that for fixed F-"* there are at most 0{y/Ym) faces F™ such 
that {F^ + fiFj") n bd (F„) / 0. For two faces F^ and FJ" we have 
R{Ff^ + (jF™) < R{Ff^) + i?(F™) < ct/VA^. Proceeding as before we see 
that there can be at most 0(1) faces of Pm which intersect (F™ + aPp). 
Altogether we have found that 

#0 = O{f^^) and #G{FP) = 0{^). (3.12) 
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Next we ask for the efficiency of Algorithm 3.1 to determine these sets. 

By the above estimate 



Let 5 -u™'* 

for the circumradius of two facets FJ" , 



be a minimal box containing F" 



we have + 'fc 



U 



■mj I 



^ I — 0{l/yff^) for each coordinate k. Again using the area bound (|3.9|) 
we see that there are at most 0(1) facets FJ" with u^'^ > Z™'* + (or 



U 



1 y 



and < u 



"1 



+ (or 



m,s ^ m,i 



"1 



n 



(cf. 



Hence the number P{P) of Lemma 3^ is of order 0(1). Next we note that 
on account of our assumptions the maximal number i^iP) of neighbours of 
a given facet is constant, if fm is large enough. Moreover, by ( |3.11 ) we 
also see that the number ry(P) is of order vTm therefore Algorithm 3.1 
determines the set G as well as the sets G{Ff^) in time 0{f^) (cf. Lemma 

O). 

Since we have already proven that #^(F™) = 0{^/f^) and that for two 
given facets FJ^, there are only 0(1) many facets FJ^ with {Fl^,Fl^, FJ^) 
e g, it remains to 'show that #{G{F{^) n G{Fl^)) = 0(1) for F/^ G ) 



(cf. Algorithm |3j). Let G ^(F™). Then we have (cf. (|3|)) 



C 



y ebd(P, 



<(p(xl),p(x2))<-- + 



/n 



2 



C5 



<(/5(y),p(x*))< 



and giF[^)ng{Fl^) 
a C5 



+ 



1,2 



The radial projection of the latter set is a spherical parallelogram with "edge 
length" 0(\l \ff^. Hence its spherical area is 0{l/fm) and together with 
(PD this shows #{G{Fp n G{F]^)) = 0(1). 

□ 

Remark. Obviously, in the cases II, HI. and IV. the sets Q{Fi), Q coincide. 
Furthermore, since (Fj — Fj) n F^ = (Fj + (—Fj)) n F^ and —Fj is a facet of 
P it suffices to determine ^(Fj), ^ only for one case. 

4 Determination of the local minima 

In this section we concentrate on the steps SI., S2. and S3, of Algorithm 



3.2. To this end let F G /Cq be a centrally symmetric polytope with n facets 



given by the inequalities 

P = {xeR^ : {ayx <bi,l<i<n}, (4.1) 

where a' G and bi G M>o. Let Hi = {x e : {a'^x = h} be the 
supporting hyperplane of the facet Fj . Since the cases I - IV of the algorithm 
can be treated " more or less" in the same way, in the following we shall focus 
only on case I and at the end of the section we shall discuss the necessary 
changes and adoptions for the other cases. 
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For six hyperplanes Hi., 1 < j < 6, the set SH,^,...,Hi^ (see step SI. and 
^) is given by all 



p3x3 



satisfying 



/ (a^^T 





\ 











(a^2)T 


{a'^y 




(a^3)T 
-(a^*)T 




K 
hi 


-(a^^)T 





(a^5)T 


K 




-(a*6)T 


) 







(4.2) 



We denote the matrix on the left ^ii,...,ig G 



Cn,...,*6, ^L-.i,^ 1 < i < 9 - rank(^,,, 



W £ 



p3x3 . 



W = Ci,. 



9— rank(ylij 



^6x9 g^j^j with suitable matrices 
we may write 



(4.3) 



We remark that the matrices Cjj....^jg, M^^^ can easily be determined by 
any program for solving systems of linear equations. 

Lemma 4.1. // rank(74j^^...^jg) < 6 then the function fHi^,...,Hi^iW) {cf. 
(2^)) has no local minimum W £ SHi^....,Hi^ with fHi-^,...,Hi^(W) > 0. 

Proof. Since the vectors a^^ correspond to facets defining hyperplanes the 
vectors (a*i , 0, 0)''', (0, a*^, O)""", (a*^, — a*'*, 0)''' G MP are linearly independent. 
Otherwise we can assume that a'^^ = a*^ = iba*^ and we get {a'''^y{w^ — 
up') = hi^ =F 7^ ^ii = ^«6- Hence the vectors {a*^, a*"', a*^} are linearly 
dependent, because otherwise rank(Aj^_,,.^jg) = 6. In the same way we find 
that {0*^,0*^,0**^} and {a*^, a*", a**^} are linearly dependent. Thus we can 
find three nontrivial vectors G M^\{0} such that 

G lin{a^^a^^a^5}^, t>2 g lin {a^S a^^ , a^«}^, G lin {a^^, a^«}^, 

where lin denotes the orthogonal complement of lin U . Now let W = 
{w^ ^uP ,w'^) G SHi^,...,Hig with det(Vl^) / 0, and let us assume that is a 
local minimum. By the choice of f * we have {w^ + Xv^ , uP + ^v^, uP + uv^) G 
SHi^,...,Hi^ for all A,/i, G M and therefore let 

g{X, fi, u) = det{uP + \v^,uP + ^jlv^ , uP + i^w^) = det(VF) 

+ A dQl{v^ ,uP ,uP) + /X det(u;^i;^,it;^) + det(it;\ u;^, v^) 
+ A/i det(t;\i;^,tt;^) + /ii/ det(^«^ f ^, w^) + Ai/ det(v\ ■w^, v^) 
+ A/ii^det(i;\t>^,'(;^). 

It is easy to see that this function has a local extremum at (0, 0, 0) if and 
only if it is constant, i.e., g{X,n,u) = det(Ty). In particular we have 
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det{v^ , w'^ , w"^) = det{'w^ , v'^ , w'^) = det(t>^, w^, w^) = 0. Since 'w^,w'^,'w^ 
are linearly independent this implies that or v'^ belongs to linji;;^}, and 
in the same way we find that v'^ or lies in linj-u;^} and or be- 
longs to linjit;^}. However, since E M^\{0} this yields the contradiction 
\det{v\v'^,v^)\ = |det(Vr)|. 

□ 

Since we are only interested in local minima W of the functions fHi^,...,Hi^ 
with fHi^,...,Hi^{W) > 0) the following we assume that rank(^j^_...^je) = 6. 
Instead of searching for the local minima of fHi-^,...,Hi^iW) = I det(VF)| it is 
more practical to look for all local extrema of the function det(l^), which 
can be parameterized by pii,...,ie : ^ M given by (see (|4.3| )) 

Pi,„„,i,(x, y, z) = det (C,,,...,i, + X ■ M^^^^ + y ■ M^^^^ + z ■ M^^,^) . 

(4.4) 

pii,...,ie is a polynomial in 3 variables {x,y,z), where each monomial has 
total degree 3 at most. In what follows we are mainly interested in general 
properties of such a polynomial, and therefore we shall write p instead of 
Pii,...,i6- So p can be written as 

P{x, y,z) = ^ aij- fc • x'y^z'', 

0<i,j,k<3, i+j+k<3 

for some scalars ,j,fc G M. The canonical first step in order to find the local 
extrema is to calculate the set V(Vp) where the gradient Vp vanishes, i.e., 

V(Vp) = {{x,y,z) eM.^ -.Vpix, y,z) =0} . 

Thus we are interested in the common roots of the partial derivatives 

pi = dp/dx = xix'^ + h{y,z) ■ X + qi(y,z), 

p2 = dp/dy = X2x'^ + hiy, z) ■x + q2(y, z), (4.5) 

P3 = dp/dz = X3X^ + l3(y, z) ■x + q^iy, z), 

where Xi G M, = Xi,2y + Xi,3Z + Xi,o, Xi,j G 1^, and qi = Xi,2,oy^ + Xi,i,iyz + 
Xifl,2z'^ + Xi,ifly + Xifl,iz + Xi,o,o, Xi,j,k G M, 1 < i < 3. As p is a polynomial 
of total degree of at most 3, V(Vp) has the following nice property. 

Lemma 4.2. Let m G V(Vp) he a local extremum of the function p and let 
C C V(Vp) be a (path-) connected component containing m. Then aff (C) C 
V(Vp). 

Proof. Let n G C, m ^ n. Since there exists a path from m to n in C 
and since the gradient vanishes on C we have p(m) = p{n). For t G M let 
s{t) = m + t ■ {n — m). The function p(s(t)) is an univariate polynomial in 
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t of degree at most three with p(s(0)) = p(s(l)) and the derivative vanishes 
at and 1. Thus p(s(t)) = p{Tn) for all t G M. Since m is assumed to be a 
local extremum there exists a t G (0, 1) such that s{t) is a local extremum, 
too. Hence we have found three points on the line s(t) where the gradient 
of p vanishes. Since all partial derivatives are polynomials of total degree at 
most two we have shown that Vp(s(t)) = for all t E M. □ 

The last lemma tells us that in order to locate all possible local extrema 
of p it suffices to find all isolated affine subspaces of Vp = 0, i.e., 

Vaff (Vp) = \ C C V(Vp) : C = aff (C) and there exists no connected 

^ (4.6) 
component U C V(Vp) with C '^U^. 

To our surprise we have not found any efficient algorithm for the deter- 
mination of the set V or V^s or, in general, for the determination of the 
local extrema of the polynomial p. Therefore we have developed an ad hoc 
method for our purposes which is based on resultants of polynomials and 
Lemma 4.2, For a detailed treatment of resultants we refer to | CL092| ] and 



[|CL098| . Here we just collect some of their basic properties. 

To this end we denote for a polynomial f G M[xi,...,Xd] by deg(f) 
its total degree and by deg(f, Xj) we denote the degree of the variable Xi 
if we consider f as a polynomial over M[xi, . . . , Xj-i, Xj+i, . . . , x^]. For a 
polynomial or a system of polynomials I C ]R[xi, . . . , x^] we denote by 
V(J) = {x £ W : f(x) = for all f G I}. Furthermore, let V^s{I) be 



the set of all isolated affine subspaces of V(X) in the sense of (4.6) and let 

Vfff (T) = {Ce Vaff (X) : dim(C) = j} , j = 0,...,d. 

The elements of V^g will be called isolated roots. Since in general is seems 
to be a hard problem to determine exactly the set Vafi (I) we only look for 
an approximation of this set. This means we want to determine a set 

Vaff (2^) consisting of finitely many affine subspace C C V(T) 
andVaff(X) C Vafr(T). ^^'"^^ 

Now let P E M[rEi, . . . , x^], j = 1, 2, be two polynomials and with respect 
to the coefficient ring R[x2, ■ ■ ■ , Xd] we write 

/•'(xi, . . . , Xd) = fi,j x\, 

1=0 

with fij E M[x2, . . . , Xd], nij = deg(P, xi) and let mi,m2 > 1. The resultant 
of / and g w.r.t. xi will be denoted by res (f-*^, f^, xi) E M[x2, . . . , x^] and it 
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is given by the determinant of the Sylvester matrix of and w.r.t. xi, 
i.e., 



^ fmi,l 
fmi 



ves{f\f 



Xl 



det 



1,1 



/m2,2 
fm2 



\ 



1,2 



/o,i 



v 



/mi,l 
fmi — 1,1 



/o,i 



fm.2,2 
fm2 



1,2 



/o,2 



/0,2 / 



(4.8) 

where the columns corresponding to f^, are repeated m2-times, mi-times, 
respectively. The definition can be extended to the case mi + m2 > 1 by 
setting res (f-*^, xi) = p if nij = 0. 

Lemma 4.3 (see [|CL092|] , pp. 150). 

i) Let {ui,U2, . . . , Ud) GM.'^ be a common root of and f^. Then res (f-*^, 

f,Xl){u2, ...,Ud) = 0. 

ii) Let u = {u2,---,Ud) € R'^^^ such that res (f^, f^, j;i)(u) = but 
fmi,i{u) or fm2,2{u) ^ 0. Then there exists a ui £ C such that 
f(ni,n) =f2(ui,u) =0. 

iii) res (f^, f^, Xl) = if and only if and f' have a common factor g E 
]R[xi, . . . , wii/i deg(g, Xl) > 1. 

In general we would like to use resultants in the following way: Let pi, 
p2 and p3 be the partial derivatives of the polynomial p (see ([4.5]) ). First 
we compute the two resultants res 1^2 = res (pi, p2, x) G M[y,z], res2,3 = 
res (p2, p3, x) G M[y, z] and then the resultant res 1^2,3 = res (res 1^21 res 2,3, y) 
€ ]R[z]. Next we determine the real roots of res 1^2, 3, and for each root z we 
determine the common roots of res 1^2(2/, 2) = and res 2,3(2/,^) = 0. So we 
get a couple of common roots of res 1^2 arid res 2,3- Again we put each pair of 
those common roots into Pi,p2)p3 and solve the three polynomials w.r.t. x. 

Now, as we shall see, if res 1^2,3 7^ (and thus all resultants are non 
trivial), then we can compute V(Vp) by this method. However, in general 
some of the resultants may vanish and hence we have to find common factors 
of the polynomials (cf. Lemma iii))- According to Lemma 4.2 we are 
mainly interested in factors corresponding to affine subspaces. Therefore we 
call a polynomial [ G M[xi, . . . , x^] linear if deg([) = 1, i.e., 

d 



i = li 



k G 
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Remark 4.1. Let f = Yh=o fi ^i^ = h + Yh=i k ■ Xi e R[xi, . . . , rr^] with 
fi G M[x2, • • • , Xd], /j G M and fm, h 7^ 0. Then I is a factor of f if and only 
if fiihxi - Vi/li,X2, ...,Xd)=0. 

Proof. Obviously, if [ is factor of f then the statement holds. W.l.o.g. let 
^1 = 1 and [ = I — xi and let q = X]i=o ^1 polynomial whose 

coefficients are recursively defined by 

Qm-i = fm and Qi = fi+i - qi+J, i = m - 2. . . . ,0. (4.9) 

Multiplication of [ and q yields l-q = f— /o+(-9o- Since /o, Iqo G M[x2, . . . ,Xd], 
[(-I, X2, - ■ ■ ,Xd) = f(-T, X2, - ■ ■ ,Xd) = we must have -/o + T • go = 0. □ 

Hence the (d— l)-dimensional affine subspaces, where a polynomial vanishes 
are given by the linear factors, and visa versa. Therefore for a set 2 of 
polynomials in M[xi, . . . , x^] let 

= {C C V(J) : C = aff (C) and dim(C) = d - 1}. 



On account of Remark |4.1| we can easily determine all linear factors and 
thus VV"^^^(f) of a given polynomial by the following procedure. 

Algorithm 4.1 (Determining a linear factor of a polynomial f). 

• For each variable xj do (without loss of generality let j = 1) 

FI. Write f as f = ^TJo fi A with fi G R[x2, ...,Xd]. 
If mi ^ do 

F2. Find d affinely independent points v? , . . . , u'^^^ G M*^"^ such that 
f(xi,uO / 0, 2 < i < d + 1. 

F3. For 2 < i < d + 1 determine the real roots Zi of the univariate 
polynomials ^{xi,u^) = 0. 

F4. For each {z2, ■ ■ ■ , Zd+i) £ Z2 x ■ ■ ■ x Zd-\-i do 

F5. Determine the solution I = {lo,l2, ■ ■ ■ ,ldy G M*^ of the linear 
system (1, n*""")/ = z-i, 2 < i < d + 1. 

F6. Let [ = -lo + xi - Ylt=2 k ■ Xi. If f(xi - i,X2, ■ ■ ■ ,Xd) = 
then { is a linear factor and the remainder f = f/[ is given by the 
polynomial q = X^I^o Qi ^1 defined in ( ^ ) . 

• If for all solutions I = {lo,l2, ■ ■ ■ JdV £ '^'^ the corresponding 
linear polynomial I is not a factor of^ then f posseses no linear factors. 



Lemma 4.4. Let f G M[xi, . . . , x^], f 7^ 0. Then Algorithm finds a linear 
factor I G M[xi, . . . ,Xd] of f and determines the polynomials f/[, if f has a 
linear factor at all. 
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Proof. Let I = Iq + Yli=i ^i^i be a linear factor of f and let us assume 
w.l.o.g. that li = —1. Let [ = /q + ^23^2 + • • • + IdXd G M[x2, . . . ,Xd\ and 
I = (^0) ^2) • • • ) ^d)^- By Remark we have f([, a;2) • • • j 3;^) = and for the 
points It* we get = f(T(M*),u*) = f((l, u*) and hence e ^i, 

i = 2,...,d + 1. By the choice of the vectors each linear system has an 
unique solution. 

Finally, we remark that the {d — l)-linearly independent points can 
be found quite easily. If the leading coefficient of / with respect to xi is 
constant, then we set = e^~^ G M'^^-'^, 2 < i < d, and u'^'^^ = 0, where 
denotes the i-th unit vector. Otherwise we perturb these points a bit. □ 

Since each polynomial can be written as an unique product of irreducible 
polynomials we can apply the above algorithm iteratively to the computed 
remainders in order to determine all linear factors of a given polynomial. 



Obviously, we can also use Algorithm 4.1 to find a common linear factor of 



two or more polynomials. Hence we have 

Corollary 4.1. Let X he as set of polynomials fj G M.[xi, . . . , x^], i £ I. 
Then there exists an algorithm which computes all common linear factors of 
the polynomials fi, i € I, and the set W^~^{Z). 

For special polynomials we have a simple test to decide whether one 
polynomial is a factor of another one. 

Remark 4.2. Let f = ZT=o fi ' ^T~' ^ ^^i^ ...,Xd],Q = YTi=o 9i ■ ^i"' e 
M[xi, . . . , Xd], rn > n > 1, fi, gi G M.[x2, . . . ,Xd] and go E M\{0}. Then there 
exists an algorithm which decides whether q is a factor of f and determines 
the polynomial f) = f/g if q is a factor of f. 

Proof. W.l.o.g. let gQ = 1. By definition g is a factor of f if and only if there 
exists a polynomial f) = X^^q" ' hi £ M.[x2, . . . ,Xd\, such that 

• f) = f. Hence by comparing the coefficients we get 

min{fc,n} 

fk= ^ gj-hk-j, k = 0,....m. 

j=max{0,n— m+fc} 

Since go = 1 the first m — n identities determine uniquely the coefficients hj , 
j = 0, . . . ,m — n, and the remaining identities can be used as a verification 
whether g is really a factor of f, i.e, whether f) = f/g. □ 

Next we describe algorithms how we can find all the isolated affine sub- 
spaces for some very special polynomials in two variables. 

Lemma 4.5. Let f = fo ■ + fi ■ x + f2 G K[a;,y] with fi G M[y] and 
deg(/j) < i. Then there exists an algorithm which computes VafT(f). 
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Proof. Since V(f) is a conic section, the set of all isolated affine subspaces 
of V(f) consists either of one or of two lines or of an isolated root or it is 
empty. Using well know formula for conic sections one can determine all of 
them. 

□ 



Lemma 4.6. Let f G M[x, y] with deg(f) < 4 and f 7^ 0. Then there exists 



an algorithm which computes a set Vaff (f) in the sense of (4.7). 



Proof. Let f = ^^^=0 /« ' ^^^"N fi ^ ^[u] and deg(/j) < i. On account of 
Corollary f4.1| we determine all linear factors of f and the set W^(f) containing 
V^g (f). Hence it remains to determine the isolated roots, i.e., (f). Since 
none of the isolated roots lies in a 1-dimensional subspace of yV^(f), we 
can divide f by the linear factors corresponding to the elements of W^(f). 
Therefore, in the following we assume that f has no linear factors. If f 
is an univariate polynomial in x, say, having a real root x, then —x + x 
would be a linear factor of f. Thus if f is an univariate polynomial we set 
Vaff (f) = WHf). Hence let deg(f,x) > deg(f,2/) > L 

If {x, y) is an isolated root then Vf (x, y) = and thus we can look for 
all isolated roots of the set {{x,y) G : f(2;,y) = fa;(x,y) = 0}, where 
denotes the partial derivative of f with respect to x, i.e., fa; = X]i=o " 
with fi G deg(/j) < i. Now let res = res(f, fa;,x) G M[7/] be the 

resultant of these two polynomials. If res 7^ then we determine V(res) 
and afterwards Wq = : f = fa; = 0, y G V(res)}. Observe, that 

or any fixed y G V(res ) the polynomial f (x, y) can not vanish, because 
otherwise it contains a linear factor (see Remark ^]^). In this case we set 

Vaff(f) = WoUWHf). 

It remains to consider the case res = and therefore we may assume 
deg(f, x) > 2. By Lemma |4.3| we know that f and fa; have a common fac- 
tor g which has positive degree in x. Even more, it is not hard to see 
that res (f, fa;) = implies that f is divisible by some 1)2, \) G M[x,y] with 
deg(f),x) > 1. Hence, if deg(f, x) < 3 then f can be written as a product of 
polynomials f^ with deg(ffc,x) = 1 which shows that f has no isolated roots. 

So let deg(f, x) = 4. Then the common factor g of f, fa; has degree 1, 
2 or 3 w.r.t. the variable x. If deg(g,x) = 1 than it is a linear factor, 
because the leading coefficient /o is a constant. Also, if deg(g,x) = 3 than 
f/g G M[x, y] is a linear factor of f. Hence deg(g, x) = 2 and it can be written 
a-s g = Yl'i=o 9i^'^~^ with gi G deg{gi) < i. Thus fa;/g is a linear factor 
and we can determine g by the following procedure: For each linear factor 
[ of fa' let g( = fa;/[. Then we have to test whether g( is a factor of f. This 



can be done with the algorithm described in Remark 4.2. So we can assume 
that we have found the factor g of f with deg(g) = 2. Let g = f/g. Then 
the problem is reduced to the determination of the isolated roots of the two 
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polynomials g and g, which can be solved by Lemma 4.5. Hence we set 



Vafr(f) = VOff(s)UVOg(g)UWi(f). □ 

Lemma 4.7. Let f, g G K[x,y] with f, g / and deg(f) < 4, deg(g) < 3. 

Then there exists an algorithm which computes a set Vafr (f , g) in the sense 
of O- 



Proof. We proceed as in the proof of the last lemma. First we determine 
all common linear factors of these two polynomials and the set W^(f, g) 



containing V^g(f, g) (cf. Corollary 4J). Then we divide the polynomials 



by these common linear factors and it remains to find the isolated roots 
V^g (f, g). To this end we may assume that deg(f, x), deg(f, y) > 1. 

Let res = res(f,g,x) G ^[y]- If res / then we determine V(res), 
Wo = {{x,y) : f = f,, = 0, y G V(res)} and we set Vaff (f) = Wo U W^f)- 
Therefore we can assume res = 0. Since deg(g) < 3 we know that if the 
polynomial g is not irreducible then it has a linear factor I, say. Using 



Algorithm 4.1 we can determine such a linear factor [ as well as the remainder 



g = g/[. Next we split our system in two systems, namely 

: f = g = and J2 : f = ^ = 0. 

Since I is a linear factor and g and f are assumed to be free of common linear 
factors the second system can be solved by substituting one variable in f via [. 
Since deg(g) < deg(g) we can apply recursively the previous argumentation 
to the system Zi. Thus, if g is not irreducible we set Vag (f, g) = V^g (Ti) U 

VOg(J2)UWi(f,0). 

If g is irreducible, and since res = 0, the polynomial g has to be a factor 
of f . Hence the determination of (f , g) is reduced to the calculation 
of the isolated roots of g. With the algorithm of Lemma [4.6| we can find 
a set Vaflf (g) C V(g) containing Vafr (g) and we set Vafr (f, g) = Vafr (g) U 
>VHf,0)- □ 

Lemma 4.8. For j = 1, . . . , 3 let f = ■ x'^ + fi ■ x + f^ £ R[x, y, z], 

fi G M[y,z], with deg(i^) < i. Let f^ = 0, / 0, res 1,2 = res (f, f, x), 
res 2,3 = res (f^, f^, x) and let res 1^2 / or res 2,3 / 0. Then for each 
L G Vaff(f^,f^,f^) there exists a M e Vag (res i_2, res 2,3) U (res 1,2, res 2,3) 
such that L C {{x, y, z) G : {y, z) G M}. 

Proof W.l.o.g. let res 1,2 / and let L G Vaff(f\f ,f )• By L'" we denote 
the orthogonal projection of L onto the plane {{x,y,z) G : x = 0}. 
Obviously, L'^ is an affine subspace and by the definition of resultants we 
have res 1^2(2/, -z) = res 2,3 (y,^;) = for all {y,z) G L'^ . Since resi_2 7^ we 
have dim(L'^) G {0, 1}. If dim(L'^) = 1, then L'^ itself corresponds to a 
common linear factor of res 1^2 and res 2,3 and thus L'^ G M^^(res i_2, res 2,3)- 



23 



Therefore let L'^ = {{y^,z^)} and so L is either an isolated 1-dimensional 
subspace or an isolated root of V(f^,f^,f^). If {y^,z^) is contained in an 
1-dimensional subspace of l^^(res 1,2 > res 2,3) the statement is certainly true. 
Thus we may assume that l{y^,z^) 7^ for every common linear factor [ 
of res 1^2 and res 2,3 and we have to show that = {y^,z^) is an isolated 
root of Vaff (res i^2i res 2,3)- Suppose the contrary and let = {y^,z^) G 
V(res 1^2) res 2,3) such that there exists a path V C V(res 1^2, res 2,3), V = 
{m* : t G [0,1]}, connecting and u^. Then V ^ conv it^}, because 
otherwise affju*^,!!^} C V(res 1^2, res 2,3) is a 1-dimensional set containing 
ti . Hence we can assume 

/2(u<) / 0, i G (0, 1), and by Lemma ii) there 
exist a*, 6* G C such that f (a*,M*) = f (a*,M*) = = f{h^M) = f 
t G (0, 1). Since fl = and /f (m*) / we get = b'^ e M. However this 
shows that L is not an isolated affine subspace of V(f^, f^, f^). □ 

Of course the last lemma also implies 

Corollary 4.2. For j = 1,2 let P = ■ + f( ■ x + £ R[x, y, z], //' G 
R[y, z], with deg(//') < i. Let = 0, / and let res 1,2 = res (f^ f , x) 7^ 
0. Then for each L G Vafr (f"*^, f^) there exists aM £ Vafr (res 1^2) UM^^(res 1^2) 
such that L C {{x,y,z) G : {y,z) G M}. 

Proof. We set f^ = f^ and apply Lemma E!^. □ 



Using Corollary f4.1| and the previous lemmas we can make the resultant 
approach practicable. 

Theorem 4.1. Let p{x, y, z) be a polynomial with total degree at most 3. Let 
V(Vp) = {{x,y,z) G : Vp = 0}. There exists an algorithm computing a 
set Vaff (Vp) in the sense of ([4?7|) . 

Proof. First we note that after some scaling, subtractions, and renumbering 
we may assume that V(Vp) is given by the following system I of polynomial 
equations (see ( [4.5D ) 

Pi = (i-x-hqi=0, 
I: p2 = [2 • a; + q2 = 0, (4.10) 

p3 = K • + [3 • X + qs = 0, 

where k G {0,1}, U G are linear polynomials and q^ G with 

deg(qj) < 2. We may further assume that deg(pi) < deg(p2) < deg(p3). 
Depending on the number of non-trivial polynomials and the number of 
variables in I we have to distinguish several cases. Obviously, if all polyno- 
mials in 4.1c vanish then we have V(Vp) = and we have to do nothing. 

(0) I consists of one or two or three polynomials in only one variable. 

Then we can determine V(Vp) by any algorithm computing the roots 
of an univariate polynomial. 
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(1,2) X consists of one polynomial in two variables. 

W.l.o.g. let V(T) = {{x,y,z) G : q3(y,z) = 0} = M x V(q3), with 



V(g3) C M . Via the algorithm of Lemma [4.5| we can determine V(q 
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and we set Vag (X) = {M x C : C G Vas (qs)}- 

(2,2) X consists of two polynomials in two variables. 

W.l.o.g. let V(X) = {(x,y,z) e : q2(y,z) = <\^{y,z) = 0} = M x 



V(q2)l3)- Using the Algorithm of Lemma [4.7| we can determine a 
set Vafr (q2,q3) C (cf. (|]^)) and we set Vag (X) = {M x C : C G 
Vaff(q2,q3)}- 

(3,2) X consists of three polynomials in two variables. 

W.l.o.g. let V(X) = {{x,y,z) e : q^(y,^) = ^^(^y^^) = q3{y,z) = 
0} = M X V(qi, q2, qs). We may assume that both variables occur in all 
three polynomials and that the polynomials are linearly independent. 
Otherwise we can reduce this case to one of the previous ones. First, 



by Corollary 4.1 we determine all common linear factors and the set 
W^(qii q2i qs) containing V^g (qi, q2, qs) and hence we may assume that 
qi; q2; qs have no common linear factors and both variables occur in the 
polynomials. Next we compute the resultant res2,s = res(q2,q3,y) ^ 
]R[z]. If res2,s 7^ then we determine W = {{y,z) : qi = q2 = qs = 
0, z £ V(res2,3)}. Observe that for a fixed z not all three polynomials 
qi can vanish, because otherwise they have a common linear factor. In 
this case we set Vaff (X) = {M x C : C e W^qi, q2, qs) U W}. 
It remains to consider the case res 2,3 = 0. Since q2, qs are assumed 
to be linearly independent and since deg(qj) < 2 the common factor 



(cf. Lemma [4.3| iii)) has to be a linear polynomial [2,3 which can be 



determined via Algorithm 4.1. Then we consider the two systems 

Xi : qi = 0, [2,3 = 

X2 : qi = 0, q2/[2,s = 0, qs/b.s = 0. 

Both systems are free of common linear factors and since both systems 
contain linear polynomials we can easily determine V(Xi) and V(l2) 
by substitution. We set Vaff (X) = {MxC : C G >V^(qi, q2, q3)U V(Xi)U 

V(X2)}. 

• In the remaining cases we first determine the set yV^(X) and there- 
fore we may always assume that the given partial derivatives have no 
common linear factors. 

'1,3) X consists of one polynomial in three variables. 

a ) Let V(X) = {{x,y,z) ^M? : [3 • x + q3 = 0} Then we may assume 
[3 7^ and the dimension of any affine subspace of V(X) is 2. Therefore 
we set V(X) =W'^{I). 
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b ) Let V(X) = {{x,y,z) £ : f = x'^ + Is ■ x + = 0}. Let 
q = [3 • [3/4 - q3 G M[y, z]. For every {x*,y*,z*) G V(X) we have x* = 
-Isiy*, z*)/2±y/q{y*, z*) and hence q{y*, z*) > 0. Uq{y*,z*) >Othen 
we can always find a neighbourhood U of (y* , z*) such that q(y, z) >0 
for all (y, z) £ U and (x* ,y* , z*) belongs to a 2-dimensional connected 
component of V(f). Hence we may set (cf. Lemma [4. 51) 



Vaff (X) = Ucev^,. (q) {(-l3(y, z)/2, y, z) : (y, z) G C} . 

(2,3) T consists of two polynomials in three variables. 

Then we may assume w.l.o.g. that the variable x occurs in both poly- 
nomials. 

a ) Let V(T) = {{x,y,z) £ : = ■ x + qi = 0, P2 = l2 ■ 
X + q2 = 0} with 7^ 0, [2 7^ and we can assume that pi,p2 are 
linearly independent. Let res 1^2 G R[y, z] be the resultant of these 
two polynomials with respect to x. Since the polynomials are linearly 
independent, deg(pi), deg(p2) < 2, and since we have assumed that 
they have no common linear factors the resultant res 1^2 can not vanish. 
By definition res 1^2 G M[y, z] is a polynomial of total degree at most 



3. By Lemma we can find a set Vafr (res 1^2) containing Vafr (res 1^2) 



and for each C G Vag (res 1^2) we consider the system 

2c : Pi{x, y, z) = 0, p2{x, y, z) = 0, (y, z) G C. 

Since C is a 0- or 1-dimensional affine subspace, Xc is a system of at 
most two polynomials in at most two variables. Hence by the previous 
cases we can determine a set Vag (Ic) C V(2c) containing Vag (Ic) and 
we set 

Vaff(2:) = W2(T) IJ Vaff(Xc). 

CeVaff (res 1,2) 



On account of Corollary 4^ we have Vas (X) C Vaff {X) ■ 

b ) Let V(I) = {(x,y,z) G ; = Irx+qi = 0, p2 = x'^ + {2-x+q2 = 
0} with [1 7^ 0. Then we can proceed as in the case above. The only 
difference is that the total degree of res 1^2 is at most 4. 

(3,3) X consists of three polynomials in three variables. 

Then we may assume w.l.o.g. that the variable x occurs in at least 
two of the polynomials, b 7^ (cf. ( |4.1C| )) and that each of the three 
polynomials contains at least two variables. Now we compute res 1^2 = 
res(pi,p2,a;) and res 2,3 = res (p2, P3, 2;)- 
a) res 1^2 = 0. 



By Lemma 4.3 we know that pi and p2 have a common factor with 



positive degree in x. Thus pi, p2 are either linearly dependent or they 



26 



have a common linear factor. If they are linearly dependent we can 
proceed as in case (2,3) b). So let [ be a common linear factor and let 
Pi = Pi/^) P2 = p2/t- Observe, pi,p2 are linear polynomials. Next we 
consider the two systems 

Ji : I = 0, P3 = 0, 

: Pi = 0, p2 =0, p3 = 0. 

Since Vafi (I) C Vafi (2^i)UVafr {I2) it suffices to consider Vafi (Ij). Since 
both systems contain linear polynomials we can reduce them to sys- 
tems in at most 2 variables which can be handled by one of the methods 
described in one of the previous cases. 

b) res 2,3 = 0. It is not hard to see that also in this case we can split 
our system in two systems containing linear polynomials and we can 
proceed as before. 

c) res 1^2 / and res 2,3 / 0. Since deg(resi^2) < 3 and deg(res2,3) < 4 



we can use the algorithm of Lemma 4^ in order to determine a set 
VafT(resi,2,i'es2,3) containing Vag (res 1,2, res 2,3). Now, for each C G 
Vaff (res 1,2, res 2,3) let 

Ic : pi{x,y,z) =0, i = l,...,3, {y,z) eC. 

Since C is a 0- or 1-dimensional affine subspace Ic is a system of 
at most three polynomials in at most two variables. Hence by the 
previous cases we can determine a set Vafr(2<:) containing Vafr(2c) 
and we set 

Vaff(X) = W2(T) J Vdc). 

CeVaff- (res 1,2, res 2,3) 



On account of Lemma 4^ we have V^s (!) C Vas (X) . 

□ 

Remark. Up to now we have only discussed polynomials arising in case 
I. of Algorithm |3.2| . Case II. can be treated completely similar to case I. In 
the cases III. and IV. we have seven hyperplanes determining the function 



fHi-^,...,Hi^i but it is easy to see that Lemma ^3] keeps true. Hence also in 
these cases we just have to examine the isolated affine subspaces of polyno- 
mials of total degree at most three. Indeed, with some extra effort one can 
show that in the cases III. and IV. it suffices to consider matrices (^ii,...,^^) 
with rank(Aj^^...^j>,) = 7. However, we note that in case IV. we have to re- 



place the 6-th hyperplane Hi^ by the hyperplane Hi^ given in Theorem 2.1 
The hyperplane Hi^^ can easily be constructed, if it exists at all. 



Altogether, using the notation of Algorithm 3.2 we have the following 
result. 
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Corollary 4.3. There exists an algorithm which computes a set Vhi^,...,Hi^ 
consisting of finitely many affine subspaces Ai, . . . , Ar of Shi^,...,Hi^^ such that 
fHi-^,...,Hi^ is constant on any affine suhspace Ai and each local minimum of 
fHi^,...,Hi^ is contained in one of the spaces A^. 

Observe, we do not determine the local minima ■M.Hi^,...,Hi^ of the func- 
tion fHi^,...,Hi^, but just a set of affine subspaces containing this set. It 
remains to check whether such an affine subspace Ai of the corollary con- 
tains an admissible lattice. To this end we use the criteria of Minkowski as 



given in Lemma 2.3, i.e., we are only interested in admissible lattices having 



a basis W £ Ai such that the conditions of Lemma are satisfied. 

Lemma 4.9. There exists an algorithm which finds for any affine subspace 
A £ Vhi^,...,Hi^ a basis W of an admissible lattice (in the sense of Lemma 



2.3) or asserts that no such lattice exists. 



Proof. For simplification we assume that Hi. = Hi, 1 < i < k. With 
r = dim A, r € {0, 3}, we may write 

^ = jvF G R^""^ : W = C + ^AiMij , 

for suitable matrices C, Mi G M^^^. Now for r > 1 and for A = (Ai, . . . , Ar) G 
R'' let W{X) = C + J2l=i ^iMi. If r = we set W{X) = C. According to 



Lemma E.3 and Theorem 2.1 we have to find a A G such that 



C bdP, (4.11) 

with j = 1 in case I., j = 2 in case II., and j = 3 in the cases III. and IV. 
Furthermore in the cases I. and II. we have the additional restrictions 

Case I: (-l,l,l)vp'(A), (l,-l,l)i4^(A), (l,l,-l)iy(A) ^ .^^2) 
Case II: (1, 1, l)^^^^-^ i P. 

Let us denote by w^(a)' 1 ^ ^ ^ the vectors of the test set U^^^y 
This means we have /c = 6 if j < 2, and k = 7 otherwise. Observe, by 
construction we know that all vectors u\^i^x) ™ supporting hyperplanes 
of the polytope. Therefore, if r = we can verify ( [4.11 ) and ( 4.12|) by just 



checking the facet defining inequalities of the polytope for the corresponding 
points Wyi/(A)) 6tc. Thus in the following we assume r > 1. Then we define 
a set A by (cf. notation at the beginning of this section) 

A={XGW:U^^,^ChdiP)} ^^^^^^ 
= {A G : u\y(^x) e H+ for ah Fi^ G AA(F'), l<i<k}. 
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Using standard methods from Linear Programming we can easily decide 
whether ^ = or we can find a point X* £ A. If ^ = then the affine 
subspace does not contain an admissible lattice. So let A* G A. In the cases 
III. and IV. or if W{X*) also satisfies ( p^ in the cases I. and II., we are 



done and we have found an admissible lattice. So let us assume that we are 



in case I. or II., A ^ $ and W^(A*) violates ( 4.12 ). The most simple way 
to decide whether there exists a A G ^ satisfying ( 4.12| ) is the following: In 
case I. we consider for i,j, I £ {1, ... ,n} the sets 



(4.14) 



Aij,i = {A G M'' : A G ^ and (-1, 1, l)wiX) e H- , 

(1, -1, l)w{x) £ H~, (1, 1, -l)vy(A) ^ Hf}. 
In case II we set for 1 < i < n 

A = {AgM'':Ag^ and {l,lA)w{X) ^ H-]. (4.15) 

Obviously, in the case I. (case II.) there exists an admissible lattice in the 
affine subspace A if and only if there exist i,j,l G {1, . . . , n} (i G {1, . . . , n}) 
such that Aij^i 7^ {Ai ^ 0). Again, using tools from Linear Programming 
we can either find a A* in one of sets Aij^i (Ai), and thus an admissible 
lattice W{X*), or we know that A contains no admissible lattices. □ 



Remark. Instead of considering the (n) feasibility problems in ( f4.14|) 
( (|4.15| )) one can argue as follows: First let us assume that we are in case 
I. and that W{X*) violates one of the restrictions in ( 4.12| ). Then we claim 



that we do not have to consider this case further: Of course, if each TV (A) 
violates (|4.12| ) for all A G ^ satisfying ( 4.11 ) then this is trivial. Hence 



suppose that there exists a A° G ^ satisfying ( [4.12 ) such that 1^(A°)Z^ is 



a critical lattice of P. Then there exists a G conv{A°,A*} C A such 
that one of the points of ( 4.12| ) w.r.t. VF(/u) lies in the boundary of P 



and the other points are not contained in the interior of P. By defini- 
tion, W{fi)Z^ is a critical lattice of P, too. W.l.o.g. let (1, 1, — 1)vk(m) ^ bdP 
and (1, —1, l)w(p) 1 1) 1)vk(m) ^ ^- Furthermore, let i/7 be a support- 
ing hyperplane of a facet of P containing (1, 1, —l)yi/(^^y Now let w^{pL) be 
the i-ih. column vector of W{^) and let W{^) be the matrix with columns 
(fi) — (fi) , u)^(/x) and —w'^{ii) + w^{iJ.). Then W{iJ.) is just another basis 
of the lattice VF(/x)Z^ and we get 

(l,0,0)vy(^) = {l,lA)w{^J.)^ (0, l,0)vy(^) = (0, 1, O)^(^) , 

(0,0, l)w{tj.) = (0, 1, (1, -'^,0)w{ii)= (l>Oi l)vy(^t)) 

(l,0,-l)vy(^) = (l,0,0),4/(^), (0, -1, l),4/(^)= (0,0, l)^(^), 

(1, 1, -l)w(,i) = (1) 1) 0)w(fj.)- 

Thus the test set P^^^ additional point (1, 1, —l)w{ii) is equivalent 

to the test set U^^^y Since W{ij,)'Ij^ is a critical lattice of P it follows from 
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the work of Minkowski (cf. |Min04], pp.27) that we shall find a basis of this 
lattice in case III. In case II. we can apply the argumentation and this means 
that in all cases it is sufficient to determine only one point of the set A. 

To sum it up, we finally have the following Algorithm for determining a 
densest lattice packing of a 3-dimensional polytope P G IC^ (cf. Algorithm 
P): 

Algorithm 4.2 (Densest lattice packing of a 3-polytope). 

Input: A polytope P G /C'^ given by the supporting hyperplanes Hi, 1 < i < 

m, or by its vertices. 

Output: A densest packing lattice of P. 

• Find the supporting hyperplanes Hi, I < i < n, of the facets Fi of 
the polytope Pq = {P — P) G /Cq (cf. (p. ID) and compute the lattice 
description of Pq. With respect to Pq do 

• For each of the cases I., II., III., IV. of Theorem do 

II. Compute the sets G{Fi) and Q with Algorithm |3.1| . 
— For three facets Fi^, Fi^, Fi^ satisfying Fi^ G ^(-^/i) Fi^ G 

g{Fl,)r^G{Fl,) do 

— For every choice of facets Fi., A < i < k, with {Fi-^, Fi^, Fi^), 
{Fi„Fi,^Fi,),{Fi„Fi,,Fi,Jeg do 

RO. IfSH,^,...,H,^^9id. dp)) do 

Rl. Determine Shi^,...,Hi^ (cf. (^)). 

// rank(Aj^^...^jj,) > 6 (cf. Lemma |4.l| and the remark after 
Theorem do 

R2. Determine a setVHi^,...,Hi^ consisting of finitely many affine 
subspaces Ai, . . . , Ar of Shi^,...,Hi^ such that fHi^,...,Hi^ is con- 
stant on any affine subspace Ai and each local minimum of 
fHi-^,...,Hi^ is contained in one of the spaces Ai (cf. Corollary 



R3. For each affine subspace A G Vhi^,...,Hi^ find a basis W such 
that W ■Ij'^ satisfies the criterion 1. of Lemma 2.3 in the first 



case and criterion 2. in the remaining cases, or asserts that 
such lattice does not exist (cf. Lemma |4.9|) . 



Among all calculated admissible lattices find one with minimal deter- 
minant. The corresponding lattice is a critical lattice of Pq and a 
densest packing lattice of P . 
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5 Densest lattice packings of the regular and 
Archimedean polytopes 

In this section we present densest packing lattices A* of all regular and 
Archimedean polytopes. Since a densest packing lattice depends on the 
representation of the polytope, and in order to make the results more trans- 
parent, we shall also give the coordinates of our representations of the poly- 
topes used for the algorithm. For a polytope P let / = (/o,/i,/2) be its 
/-vector, i.e., fi is the number of i-faces. Futhermore we shall use the fol- 
lowing abbrevations: Let r = (1/2) (1 + \/5) and let Ds be the fcc-lattice 
with basis (1,1, 0)t, (1,0, 1)t, (0,1, 1)t. Finally, let Phc be the so called 
rhombic triacontahedron given by 



rtc 



\TXi\ < 1, 



1 T r + 1 

— xi H — xo H x-i 

2 2 2 



r T + 1 1 

—xi H Xo H — xs 

2 2 2 



< 1, 



< 1, 
1 



1 r 

-X^ H Xo H 

2 2 2 



(5.1) 



For the identification of the polytopes we shall use the Wythoff symbols. 
• Tetrahedron, 3| 23, / = (4,6,4) (cf. Poy70| |). 



Pt = {x ■.Xi+X2 + X3< 1, 



A*{Pt) 
S*{Pt) 



Xl + X2 - X3 < 1, 



-Xl - X2 + X3 < I, 
Xl- X2- X3< l}. 



1 1 



0.367346938. 
(8,12,6). 



6' 6' 6 

18 
49 

• Cube, 3| 24, / 

Pc = {x :\xi\ <l} , A*(Pc) = 2Z^ 6*{P, 

• Octahedron, 4 | 2 3, / = (6, 12, 8) (cf. f|Min04|| ). 

Po = {xeR^ : \xi\ + \x2\ + Ixsl < 1} , 

A*(Po) = 2 



,1)T Z 



1. 



(- - -y 



- -y (-- - 

3'2^ 2'6' 



1. 



1 8 



0.947368421. 



• Dodecahedron, 3 | 2 5, / = (20, 30, 12). 

Pd = {x eR^ : \txi\ + \x2\ < 1, \tx2\ + |x3| < 1, Irxsl + < 1} , 

A*((l + T)Pd) = 2D3, 6*{Pd) = ^ 0.904508497. 
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Icosahedron, 5 | 2 3, / = (12, 30, 20). 

Pi = {xeR^: |xi| + \X2\ + \X3\ < 1, |rxi| + |(l/r)2;3| < 1, 

\tx2\ + \{l/r)xi\ < 1, Irxal + |(1/t)x2| < l}, 
A*((l+r)Pi) = 2{w^{x),w'^{x),w^ix))Z^ where 

/{-f - fV5)x^ + if + fV5)x - ^ - IV5\ 
w\x)= I (_i _ 1 ^) a; + 1 + 1 ^ 

V (f + f %/5)x2 + (-f - 8^/5)x + f + I , 

/(-i V5-f)x^ + {§V5 + f)x-l-§V5\ (5.2) 

V (-iV5-f)x2 + (|V5 + f)x-A^ 
u;^(x) = ((iV5 + 1)^-2- \/5, X, 0)\ 

and X is the unique root with x G (1, 2) of 
1086 + (-1063 - lllVb) x^ + (15\/5 + 43) x + 102 + 44\/5, 

(p.\ ^("^ ~^ ^) ^ Q §2635T445 

* I det(t(;i(x), it;2(x), i(;^(x))| 

Cuheoctahedion, 2 | 34, / = (12, 24, 14) (cf. iHoy70(| ). 

Peo = G M3 : X G n 2 • P4 , 

45 

A*(P,o) = A*(Pt), <5*(Peo) = — « 0.918367346. 

Icosidodecahedron, 2 | 3 5, / = (30, 60, 32). 
P,d = {xGM3:xGP,nP4, 

14 + 17r 

A*((l + T)Pid) = 2D3, 6*{Pu) = ~ 0.864720371. 

Rhombic Cubeoctahedron, 34| 2, / = (24,48,26). 

Prco = {x G : |xi| + |X2| < 2, |xi| + |X3| < 2, |X2| + |X3| < 2, 

and X G \/2 • Pc n (4 - V2)Po}, 
A* (Prco) = 2D3, 5* {Prco) = ~ ^ 0.875805666. 



Rhombic Icosidodecahedron, 3 5 | 2, / = (60, 120, 62). 

Prid = {x G X G (3r + 2) • P^c n (4r + 1) • Pi n (3(1 + r)) • P^}, 
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_ / (r-l_ 7 9r + 4y /9r + 4 r-l_ 7 



7 9t + 4 T - 1 



2' 4r + 2' 4r + 2 



<5*(P„d) = ^ f , ~ 0.804708487. 
36r + 15 



Truncated Cube, 23 | 4, / = (24,36, 14). 
P^^^ = |x G : X G Pc n (1 + \/2)Po} , 

A*(Pt,e) = 2 ((1, -a, Oy, (0, 1, -ay, {-a, 0, 1)^) Z^, a = (2 - ^/2)/3, 

= w 0.973747688. 

5 + 3^2 

Truncated Octahedron, 24 | 2, / = (24, 36, 14). 
A*(Pi™) = 2 ("(1, 0, 0)T, (1, 1, oy, (i i, -i)!") Z3, 5*(P,,,) = 1. 



Truncated Dodecahedron, 23 | 5, / = (60, 90, 32). 

Ptrd = {x G X G (1 + r) • Prf n ((7 + 12r)/(3 + 4t)) • PJ, 

1 5-1- _|_ If) 

A*(Pi,rf) = 2Z?3, <5*(Pt.d) = - • -—— ~ 0.897787626. 

4 6r — 3 

Truncated Icosahedron, 2 5 | 3, / = (60, 90, 32). 

Pt„ = {x G : X G (1 + r) • P n (4/3 + r) • Pa}, 
A*(Pt„) =A*((1 + T)P.), 

''^^^"^= 10S.\,em^^ ^ 0.7849877759. (cf. (H 
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Truncated Cubeoctahedron, 2 34|, / = (48,72,26). 

Ptrco = (x G : laril + \x2\ < 2 + 3^2, \x2\ + \x3\ < 2 + 3^2, 

k2| + |x3| < 2 + 3^2 and (2^2 + l)Pc n (3^2 + 3)Po}, 

A* (Ptrco) = 2 ((2^2 + 1, -2^2 - ^ + a, 2^/2 + ^ - a)T, 

where a = -V33(V2 + 1), 
6 

99 231 ^ 2835 ^ 6615 

S (Ptrco) = * \/66 * V33 + *V2 0.849373252. 

V trcoj ^ 1984 ^ 992 ^ 1984 

Truncated Icosidodecahcdron, 2 35 |, / = (120, 180,62). 

Ptrid = {xeR^ :xe(5T + A)- Prtc n (6r + 3) • n (5(1 + r)) • Pa] 

^*(lPtrid) = ^Ds, S*(Ptrid) = 1^ + ^- 0.827213595. 



Truncated Tetrahedron, 2 3 | 3, / = (12, 18, 8). 

Ptrt = {x e : X e 5 • n -3 • Pt}, 

/24 42 42 \, 

^•(^«' = K(3-^-3'^-(^--3-3>^-(-3-3'-^>T' 

207 

(5*(P4^t) = — 0.680921053. 
Silub Cube, I 2 3 4, / = (24, 60, 38) . 

Let Psc be the snub cube such that the 6 quadrangle facets he in the 
hyperplanes {x G : Xj = ±1}, 1 < i < 3, and let y* be the unique 
real solution of + + y = 1. 

^*(Psc) = 2 • ((1, 0, 0)T, (0, 0, 1)T, (i, 1-1, -1)T^ , 

<5*(P,e) = ^ + + ^(y*)' ~ 0.78769996. 

Snuh Dodecahedron, | 2 3 5, / = (60, 150, 92). 

Let Psd be the snub dodecahedron such that the 12 pentagonal facets 
lie in the supporting hyperplanes of the facets of the dodecahedron 
(l + r)Pd. 

A*(P,a) = 2Ds, 6*(Psd) = vol(P,d)/16 « 0.788640117. 
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Theorem 5.1. The above list contains the densities of a densest lattice 
packing of all regular and Archimedean polytopes. 



Proof. Algorithm 4.2. □ 
Remarks. 

a) In order to get exact values for the densities as given in the above list 
we first use a numerical implementation of Algorithm |4.2| , by which 
we determine an optimal selection Fi-^ , • • • , Fij^ of facets corresponding 
to a critical lattice. Then for this special choice of factes we carry out 
the steps Rl., R2. and R3. with the symbolic computer algebra system 
Maple V Release 5. 

b) The last three polytopes Ptrt, Psa Psd are not centrally symmetric and 
thus one has to calculate the difference bodies first. The difference 
body of Psc is a polytope with 74 facets, whereas \{Psd — Psd) has 
already 182 facets. We also have computed the densest packing lat- 
tice of the dual polytopes of the Archimedean polytopes, the so called 
Catalan polytopes. Thereby we had to determine packing lattices of 
polytopes with more than 380 facets. The CPU-time for the determi- 
nation of the densest packing lattices of all regular and Archimedean 
polytopes is about 5.5 hours on a PC with a 266 Mhz-Pentium II pro- 
cessor. 

Acknowledgements. We would like to thank Utz-Uwe Haus, Achill Schiir- 
mann and Chuanming Zong for helpful comments. 
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